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ABSTRACT 

I calculate the evolution of Be and B abundances produced by cosmic rays generated by massive 
stars in the pregalactic phase of the universe. The inputs for calculation, i.e. the star formation 
rate and the nuclear abundances of cosmic rays, which I assume to be the same as those of the ISM, 
are taken from the results of a detailed cosmic chemical evolution model with its parameters best 
fitted from several items of observational information including an early reionization of the IGM by 
z ~ 15. I found that when the ^Li plateau abundance observed in metal-poor halo stars originated in 
the pregalactic cosmological cosmic ray nucleosynthesis, Be and B simultaneously produced with ^Li 
amount to the lowest levels ever detected in metal-poor halo stars. It is desirable to observe Be and 
B abundances in metal-poor halo stars with [Fe/H]< —3 in order to elucidate the possibility of early 
^LiBcB production by pregalactic supernova cosmic ray nucleosynthesis. 

Subject headings: cosmic rays — cosmology: theory — nuclear reactions, nucleosynthesis, abundances 
— stars: abundances — stars: Population II — supernovae: general 



1. INTRODUCTION 



The lithium abundances 

halo stars (MPHSs) show 

tion of metalli city ( ISoite fc Spiti 119821: iRvan et all 



observed in metal-poor 
a plateau as a func- 



tion or metaiiicity i i&pite &: opita iKvan et ai.l 

20001: iMelendez & Ra mirezj |2004 lAsplund et al.l l200a 
Bonifacio et'aiTl2007c iShi et al.ll2007f ) at ^LI/hT 
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xlO"^". The prediction by the standard big bang nu- 
cleosynthesis (BBN) model of ^Li abundance which is 
the main lithium isotope observed in MPHSs, how- 
ever, indicates a factor of 2 — 4 larger value, when the 
baryon-to-photon ratio deduced from parameter fits to 
the temperature fluctuations of cosmic microwave back- 
ground (CMB) radiation measured with Wilkinson Mi- 
crowave Anisotropy Probe fW MAPl (ISoergel et al.l l200a 
l2007f) is used. For example, ICoc et al.l ()2004[ ) derived 



^Li/H=(4.15+o:45) X 10-1° with the baryon-to-photon ra- 
tio 77 = (6.14 ± 0.25) X 10-1°. This discrepancy between 
the observations and the BBN-t-CMB prediction of ^Li 
abundance is a problem, which indicates some destruc- 
tion process of ^Li. Recently, a c omple x but consistent 
theory is suggested bv lPiau et al .^ (2006) . In their theory, 
an extremely high efficiency of engulfment of baryons in 
a first generation of stars results in a half or one third of 
destructions of D and Li isotopes. Population II (Pop II) 
stars are made of a mixture of ejecta of supernovae (SNe) 
of the first stars and unprocessed material of the primor- 
dial composition, and experience a depletion of lithium 
isotopes in their atmospheres, which are observed. There 
are two important conditions. The ejecta of SNe of the 
first stars needs to mix with the pure BBN-composition 
matter at 2.5 < [Fe/H]^ in order to form the lithium 
plateau. An infall of the intergalactic medium (IGM) 
needs to proceed after the formation of Pop II stars in 
order to be consistent with observations of deuterium 

^ Research Fellow of the Japan Society for the Promotion of 
Science. 

2 [A/H]=log(A/H)-log(A/H)0 



abundance. 

Recent spectroscopic observations of MPHSs also pro- 
vide abundances of ^Li isotope. They indicate a likely 
primordial plateau abundance, similar to the well known 
^Li plateau, of ^U/R=6 x IQ-^^^ which is about 1000 
times as large as the BBN prediction. Since the standard 
Galactic cosmic ray (CR) nucleosynthesis models pre- 
dict negligible amounts of ^Li abundance with respect to 
the o bserved plateau level at [Fe/H] < —2 (e.g. iPrantzosI 
I2006D . this plateau causes another problem, which indi- 
cates some production process of °Li. 

Several candidates for early °Li production mech- 
anisms have been suggested. The non-thermal nu- 
clear reactions triggered by the decay of long-lived 
particles i s one poss i bility o f t he non-standard pro- 
cess ( JedamzikI |2000[ l2004allbl. 120061: iKawasaki et all 
20051 iKusakabe et al.l [20061 : ICumberbatch et al.l 120071) . 



Pospelovl (|20'06i r suggested the exotic nuclear reaction of 



^iiex{d,X~)^Li to make abundant ^Li, where X~ is a 
negatively charged massive particle assumed to decay in 
the early univ erse, and Hex i s the state that has "^He 
bound to X~. ISuzuki fc Inoud ()2002D suggested an a -I- a 
fusion reaction with a particles accelerated by hierar- 
chical structure formation shocks, thought to have been 
operative at the Gala xy formation epoc h. 

As a possibility, Rollinde et al.l ()2005[ ) have calculated 
the ^Li production by an initial burst of cosmological 
cosmic rays (CCRs) to show that this process through 
a + a fusion can acc ount for the ^Li plat eau without 
overproduction of ^Li. iRoUinde et al.l ()2006l ) applied the 
CCR nucleosynthesis to a well grounded detailed model. 
They derived the total CCR energy as a function of red- 
shift from a star formation rates (SFRs) in mo dels of 
cosmic chemical evolution of iDaigne et al.l (120061 ). which 
are made to reproduce the observed cosmic SFR, SN II 
rate, the present fraction of baryons in structures, that 
in stars, the evolution of the metal content in the inter- 
stellar medium (ISM) and IGM, and early reionization 
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of the IGM. As a result they found that the pregalac- 
tic production of the ^Li in the IGM via Population III 
(Pop III) stars can account for the ^Li plateau without 
overproduction of ^Li. 

The CR nucleosynthesis also produces Be and B by 
spallation reactions between CNO nuclei and p and a 
particles. If this OCR nucleosynthesis scenario of ^Li 
production leads to overproduction of Be and B nuclides 
against observations in metal-poor stars, it cannot be 
achieved in the real universe. This study is devoted to 
checking if the CCR nucleosynthesis as a mechanism of 
^Li production is consistent with observations of Be and 
B in metal-poor stars, and finding constraints on condi- 
tions of the CCR nucleosynthesis. 

Abundances of Be and B are observed in MPHSs. 
A trend of ^Be abundance is found such that Be in- 
creases lin early as Fe during th e course of Galactic 
evolution ([Boesgaard et al.l Il999f) . The first report 
on observation of beryllium in two very metal-poor 
stars with the Ultraviolet and Visible Echelle Spectro- 
graph (UVES) mounted on the ESQ VLT Kueyen tele- 
scope (jPrimas et al.ll2000a[) said that the trend of beryl- 
lium with metallicity keeps decreasing at lower metal- 
licities with no evidence of flattening. iPrimas et al.l 
()2000br ) found that the very metal deficient star G 64- 
12 ([Fe/H]=-3.3) has Be of log (Be/H)= -13.10 ± 0.15 
dex, which is significantly higher than expected from the 
previous trend, and claimed that this high [Be/Fe] ratio 
may suggest a flattening in the beryllium evolutionary 
trend at the lowest metallicity end or the presence of 
dispersion at ea rly epochs of the Galactic ev olution. On 
the other hand, iBoesgaard fc Novickil (|2006f) found that 
G64-37 with [Fe/H]=— 3.2 has a Be abundance which is 
consistent with the Be-Fe trend, and suggested that dif- 
ferent Be values are indicative of a Be dispersion even at 
the lowest metallicities. 

B abundances in metal-poor stars have been estimated 
by observations with the Goddard High Resolution Spec- 
trograph (^llRS}_on_boa£d th£_^^ 
(HST) ([Du ncan et al.' 199?: 'Garci a Lopez et all 119981 : 
IPrimas eTa l. 1999; Cunha et al. 2000-). The boron abun- 
dances are also found to show a linear increase with a 
slope of ~ 1 with respect to metallicity, and there is no 
signature of a primordial plateau abundance. 

These trends of Be and B abundances are ex- 
plained by Galactic CR nucleosynthesis models of dif- 
ferent types. One is the acceleration of metal-rich 
CRs, probably freshly synthesized matter at SNe fol- 
lowed by the primary reactions between CR acceler- 
ated CNO nuclides and interstellar nuclides p and a, i.e . 
[CNO]cR+[HiSM ^[LiBeB] (e.g. iRamatv et al.|[T997l ). 
This m echanism leads to the same-rate increase of BeB 
and O. [Fields et al.l ()200Q) have suggested that the sec- 
ondary reactions between CR accelerated pa and inter- 
stellar CNO i.e. [pa]cR+[CNO]isM ^[LiBeB], perhaps 
without any contribution of the primary reactions, would 
explain the BeB to Fe trend if the ratio O /Fe increases to- 
ward low metallicity. iValle et "aTl (|2002f ) have shown that 
the multi-zone (halo, thick disk, and thin disk) Galactic 
evolution model including only the secondary reactions 
can reproduce the linear trend without fine-tuning. 

In this work, I adopt e d Mo del 1 and the rapid burst 
model of iDaigne et al.l (|2006( ) for the cosmic chemical 
evolution model. This chemical evolution model and 



other inputs, as well as calculation of light element evo- 
lution are explained in Sec. 2. I present results of the 
CCR nucleosynthesis in Sec. 3, and discuss this study in 
Sec. 4. I summarize the CCR production of LiBeB in 
Sec. 5. 

2. MODEL 

2.1. Cosmic SN Rate 

I adopt the cosmic SFRs in a chemical evolution model 
given by Daig ne et al.l ()2006D . I take Model 1 and the 
rapid burst model, which include formation of stars with 
masses between 40 and 100 Mq in the early phase of the 
universe. As their best model, a parameter, the min- 
imum mass Mniin of dark matter halos of star-forming 
structures is determined to be 10^ Mq for chemical evo- 
lu tion of structur e s. Th e adopted models are the same 
as iRoUinde et all (|2006D use. The birthrate function is 
given by 

B{m,t,Z) ^ Mm)Mt) + Mm)MZ), (1) 
where m,t, Z are the mass of star, the age of the universe, 
and the metallicity, and 4>2 are initial mass functions 
(IMFs) of the normal and massive-only component of 
stars, respectively, Vi E^^d ^2 are SFRs for respective 
components. The normal component is given such that 
its mass range is from 0.1 to 100 Mq. The massive com- 
ponent is active only at high redshift, and its mass range 
is from 40 to 100 Mq. The IMF of both modes is given 
by a power law of mass with an index like the Salpeter's, 

0,(m) cxm-(i+^), (2) 
with X — 1.3. The amplitudes of the IMFs are normalized 
respectively as 

/ dmm(j)i{m) = 1, (3) 

where m-inf and msup are the lower and upper ends for 
the mass range of each mode. 
The normal-mode SFR is given by 

7/ii(i) = Z/iMstructexp(-i/Tl), (4) 

where vi = 0.2 Gyr ~^ describes the efficiency of the star 
formation, and Mstruct and n = 2.8 Gyr are mass of the 
structure and timescale, respectively. 
In Model 1, the massive-mode SFR is given by 

i>2{t) = I/2MiSMexp(-ZiGM/2'crit), (5) 

where 1^2 = 80 Gyr ~^ is related to the efhciency of 
star formation. AfisM is the baryonic mass of the ISM. 
ZiGM is the metallicity of the medium between the col- 
lapsed structures (identified as the IGM), and Zcrit = 
10~^Zq determines the effective epoch of the end of 
Pop III star formation. In contrast, in the rapid burst 
model, the massive mode star formation occurs as an in- 
stantaneous event at redshift z — \Q. Since the SFR in 
the model has an instantaneous bump and I f o und d if- 
ficulty reading it from Fig. 13 of iDaigne et al.l (|2006f ). I 
fixed the SFR so t hat it is ^ 2 0M(7)yr~^Mpc~^ during 
3 X 10^ yr (R,ollind e~erani2006( l. 

2.2. LiBeB Production in the Homogeneous Universe 

I calculate abundances of light elements (LiBeB) pro- 
duced in the homogeneous universe through the interac- 
tion between fast nuclei accelerated in the SN shocks and 
background nuclei. The picture of CCRs, i.e., the total 
kinetic energy of CRs and th e propagation is identical to 
that of iRollinde eTall (|2006f ). 
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2.2.1. CCR Energy and Its Spectrum 

The total kinetic energy given to accelerated CRs by 
SN explosions is 

esn{z) = {i + zf 

X / dm'^(j)i{m)^.i{t ~ T{m))£cR{m), 

Jmax(8M0,md(t)) 

(6) 

where m^it) is the mass of stars with lifetime t, T{m) 
is the lifetime of a star of mass m, and i?cr(?ti) is the 
ener gy imparted to CRs per S N with its initial mass 
m. iDaigne et all (|2004 120061 ) adop t stellar hfetimes 
Tim) from iMaeder fc Mevnetj (|1989l ) f or interniediate - 
mass stars (< 8 Mq), and those from iSchaererl ()2002f ) 
f or massive stars (8_Mq < m < 100 Mq). 

iRollinde et all (|2006[ ) give Scnim) after some assump- 
tions to calculate the total kinetic energy EsNiz). Core 
collapse SNe are supplied with energy by the gravita- 
tional collapse of cores. Almost all of the energy gener- 
ated from core collapse, i?cc is transferred out by neu- 
trinos. Only 1 % will be give n to the energy of SN ex- 
plosions. IRollinde et all (|2006D use a parameterization 



£cR{m) 



100 



(7) 



where e is the fraction of SN explosion energy imparted 
to CRs. Assumptions related to Ecc are as follows: 
Every star of mass to > 8 Mq explodes as SN. Stars 
of mass 8 Mq < to < 30 Mq make neutron stars of 
mass 1.5 Mq at core collapses and Eqc = 3 x 10^^ ergs. 
Stars of mass 30 Mq < to < 100 Mq become black 
holes with masses ver y similar to their helium core 
mass (|Heger et al. l l200l The mass of the helium core is 
Mne = 13(to - 20 Mq) (|Heger fc Wooslevll2002D . In this 
case Ecc is proportional to the mass of the black hole, 
and Ecc = 0.3Mhc is assumed. 

I calculate light elem ent production i n the f ormal ism of 
iMontmerld (flOTl and lRollinde"eraI1 (j2005l . l200l . The 
proper source function of SN CRs Qi{E,z) in a rapid 
burst at redshift Zs (corresponding time t^) is defined as 

Q,{E, z) = (1 + - ts) 

P 

[(GeV/nucleon)-i cm^^g-i], (8) 

where E and /? are the kinetic energy and velocity of CRs, 
respectively, and I define the CR injection spectrum of 
nuclide i as 

SdE.Zs) ^ KfJ^izs)—, ^- r-^, (9) 

where Kf?^ is the ratio of number abundance of i to that 
of p, i.e. i/p of CRs, and Eq — 938 GeV is the nuclear 
mass energy per nucleon. The amplitude of the source 
function is set so that SNe from both the normal mode 
(Pop II) and the massive mode (Pop III) stars supply the 
CR energy, 

£sn{z)= / Ey^Q,{E,z)dE. (10) 

I take i;n,in = 0.01 MeV, ^^ax = 10^ GeV , and CR spec- 
tral index 7 = 3 as iRolhnde et all (|2006f ) do. Since the 



evolution of CR confinement by a ma gnetic field is diffi- 
cult to estimate (jRollinde et al.ll2006l ). as a first step, I 
assume that the CR confinement is ineffective in the early 
universe, so that all CRs generated by SNe in structures 
immediately escape from structures to the IGM. In this 
case, there is uniformity of the CR density in the uni- 
verse. 

2.2.2. Primary Light Element Production by SN CRs 

I define the number density of a CR species i of energy 
E at redshift z as Ni{E,z) [in cm"'^ (GeV/nucleon)"^]. 
In order to delete the volume changing effect by cosmic 
expansion, I define the relative number abundance to 
that of the background proton nYi{z), 

N,,n{E,z)=N,{E,z)/nn{z). (11) 

The transport equation for A'i^, under the isotropic con- 
dition, is (Mo ntmerlc 1977l ) 



dt 

where h{E, z) 



(12) 



(dE/dt) 
-i] for 



is the energy loss rate 
[(GeV/nucleon) s~^] for cosmic expansion or ioniza- 
tion, and Ty){E,z) is the lifetime against destruction. 
Qi^}i{E, z) = Qi{E, z)/n-ii{z) is the normalized comov- 
ing source function. 

The expansion loss and ionization loss are ex- 
pressed in a product of energy-dependent term and a 
redshift-dependent one, b(E, z) = -B (E)f(z). These 
terms are given in IMontmerld (|1977D . The redshift- 
dependent term of the expansion loss is /e = (1 + 
z)~^\dz/dt\Hl^^ , where Hq is the Hubble constant. I 
assume the standard ACDM model with its para meters 
from WMAP three-year data (|Spergel et al.ll2007f) ^ , h = 
i7o/(100 km s"^ Mpc"^) = 0.704, fl^h'^ = 0.022, = 
0.27, Oa = 0.7 3. The ionization loss ra te is from the fit- 
ting formula in lMeneguzzi et al.l (|1971 ^ with the number 
fraction of ^He, He/H=0.08. Tq = (nH(z)crD,,/3)~^ is es- 
tima ted with the n uclear destruction cross section (Td,?; 
from' ReevesI (|1974| ). 
I define z''{E,E',z) as in lMontmerli (|1977D . 



dz* 
'dE 



1 



dz 
'di 



dz* 
dz 



(13) 



B{E)f{z) 

A physical interpretation is that a CR particle with en- 
ergy E at redshift z had an energy E'{> E) at red- 
shift z*{E,E',z) before experiencing energy loss. Thus 
z*{E,E,z) = z is satisfied. CCR particles with en- 
ergy E at z originate in those with E[.{E, z, Zs) at Zs- 
E'g{E, z, Zs) satisfies an equation, z*{E, E'^, z) = Zs- 
z*{E,E' , z) is obtained by integrating Eq. ([13]) applying 
the greater loss process to b assuming that the process 
with the greater rate of b{E, z) is dominant all the way 
from redshift z* to z. 

The transfer equation is solved (jMontmerld 11977^ 1 to 
obtain the CCR energy spectrum from a CR burst at 

Zc, 



^i.B.(E,z,Zs) = C{z^ 



exp{-aE,E'^,z)) 



dz 
It 



\b{E,Zs)\ \dz*/dE'\ 



-,(14) 



http: / /lambda. gsfc.nasa.gov 
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where $i.H(£', z, z^) = ^i{E,z,Zs)/n-iiiz) is the normal- 
ized flux of i per comoving volume with ^i(E, z, Zs) = 
(3Ni{E, z)z^, fi and (3' are the velocities corresponding to 
energy E and E'^ , respectively. n|\ is the present average 
number density of protons in the universe. ^ is an effect 
resulting when the nuclear destruction is considered, and 
given as 



£.{E,K,z)^ 



E' 



dE" 



Ie HE", z*iE, E", z))Tj,{E", z*{E, E" , z))\ ' 

(15) 

After analysis with Eq. one can estimate 

\dz*/dE'\E'=E'^ = \b{E'^,Zs)\-^\dz/dt\^^^^ and find an 
expression for $i h 

(16) 

The production rate of light element I of energy E, 
produced at redshift z is given by 



dt 



dNi^-a{E,z,z,) y-^ f / 

—=2_^ a^j^i{E,E )nj(z)^i^ii[E ,z,Zs)dE 

= Y,ja..j^i{E, E')K^^^{z)^,{E', z, z,) dE' , 



(17) 



where (Jij^i {E, E') is a cross section of a process between 
a CR nuclide i with energy per nucleon E' and a back- 
ground species j to make a given light element / with E, 
and nj{z) and K}j^^{z) are background number abun- 
dance of a nuclide j and number ratio of j to proton, 
respectively. When the destruction of the light element 
I after production is neglected, the total production rate 
is calculated as 



dNuYi{E,z,z,) 



E 



dt 



i^j^W I ^r;UE')W,z,z,)dE', (18) 



dE 



where (j\°^i{E') is the total cross section of a reaction 
i + j ^ I + X , with any X. I adopt cross sections from 
iRead fc Violal (| 19841 ) , and particularly for the a -f a re- 
action, exponential-plus-constant cr oss section for I =^L i 
and exponential one for / ^■^Li from lMercer et al.l()2001h . 
The resulting light element abundance is obtained as the 
CR production added to the BBN yield. The yield by 
CR nucleosynthesis is the integration of those produced 
at z' from CRs generated at Zg over z' and Zs, thus 

^/ IGM V^ / BBN 



dz. 



dt 



dZf, 



dz' 



dt 



dz' 



^/f]«^(z') J <jl°UE')ME',z',z,)dE'. 



(19) 



2.2.3. Secondary Light Element Production by SN CRs 



I also calculate the LiBeB production in 
the universe by the secondary process, i.e., 
[pa]cR+[CO]isM ^[LiBeBjisM- Since the C and O 
abundances of the ISM in structures are about two 
orders of magnitud e highe r than those of the IGM (see 
Fig. 11 in Daigne et al.l 12006). the secondary LiBeB 
production in the IGM is not important. I expect that 
the LiBeB abundances in the ISM are enhanced by 
a contribution of the secondary process. In fact, the 
reactions of [pa]cR+[CO]isM ^[LiBeBjisM make hght 
elements in the ISM and the mass accretion to the 
structures from the IGM dilutes the ISM abundances 
in the framework of this model involving a hierarchical 
structure formation. Note that from the assumption 
that the confinement of CRs by a magnetic field is 
ineffective, the CRs do not stay in the structures. 

The light element abundances produced by the sec- 
ondary reactions are then given with a parameter: the 
fraction of baryons at redshift z which are in structures 



/(z) = Si-" 



Pdm 



(20) 



where fps(M, z) is the distribu tion function of halos 
taken from the iSheth fc Tormenj (1999) modification to 



the Press-Schechter function (iPress fc Schechter 
converted into the mass function Jenkins et al 



1974 ) 



mm 



by a code provided by A. Jenkins (2007, private commu- 
nication). I assume that the primordial power spectral 
slope is n=l, the rms amplitude for mass density fluctua- 
tions in a sphere of radius 8 Mpc is Ug — 0-9, and the 
iBond fc Efstathioul (|1984[) fit to the transfer function for 
cold dark matter is used in generating a mass function. 
Pdm is the comoving dark matter density of the universe. 

The light elements made by the secondary process are 
contained in the structures that grow gradually. The 
abundance of a light element in the structures is then 
given by 



ISM 



(^) 
IGM 



/(^) 



dz. 



dt 



dz. 



dz' 



dt 



dz' 



x}2K]l^{z')f{z') 



3. RESULTS 



^i{E')<^,iE',z',Zs)dE'. 



(21) 



I calculate the light element production in the uniform 
universe by CCRs (i.e. neglecting an inhomogeneity of 
CCRs). I consider only processes between the acceler- 
ated C Rs with the a b undan ce patterns of the structures 
in the iDaigne et al.l (|2006l ) model and the background 
IGM abundance which I assumed to be of the primor- 
dial abundance. I set the primordial helium abundance 
Hc/H=0.08. 

I show the result of light element evolution in Model 1 
in Fig. [TJ Solid lines correspond to the case where both 
the normal and massive modes are included, and dashed 
lines correspond to the case in which only the normal 
mode is included as the energy source. When the SN 
energy is totally given to the CR acceleration (e = 1), 
6Li/H=2.0 X 10"" at z = 3 is derived. In Fig.[Tl e = 0.31 
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is assumed so that ^Li/H=6 x 10^^^ at z = 3 results that 
is the observed abundance level in MPHSs. In the case 
that the CRs are energized by only the normal mode 
star formation, the energy fraction e = 0.73 is needed to 
realize the observed ^Li abundance. The time evolution 
of light element abundances in the rapid burst model 
is shown in Fig. [21 The energy fraction given to the 
CR acceleration is assumed to be e = 0.029, and the 
abundance of ^Li gets ^U/Il=6 x 10"^^ at z = 3 These 
results are very similar to that in lRoUinde et al.l ()2006f ) 
(See their Fig. 2). In Model 1, ^Li is produced gradually 
with decreasing redshift, while in rapid burst model, it 
is immediately produced by most part at the burst of 
star formation. In my calculation, the ^Li production 
at high redshift t ends to be slightly more efficient than 
in lRollinde et al.l (jlOOfi.). This difference might be caused 
by different numerical calculation of transport equation 
for nuclides. 

Figure [3] shows the yields by CCRs generated at Zg 
per second in Model 1, i.e. A{l/H)/ Atg. The calcu- 
lated results show that ^Li and ^Li are produced mainly 
by the a + a fusion reaction, and Be and B produc- 
tion has the strongest contribution from the 0+p (and 
C+p) spallation processes. I have taken the helium abun- 
dance He/II=0.08, and the O abundance in the ISM of 
structures in Model 1 is r oughly constant fr om z ^ 
to z ~ 20 (see Fig. 11 of lDaigne et"alll2006l ). The in- 
jected energy density in CRs is si noothly increasing a s 
a function of redshift z (Fig. 2 of lRollinde et al.ll2006f ). 
This figure therefore reflects the abundance of seed nu- 
clides and injected energy density. The CRs generated 
at Zs ~ 3 have little time to react with background nu- 
clear species, so that the yields go to zero, and a decrease 
in all yields at Zs 5, 25 reflects the sh ape of CR energy 
density (Fig. 2 of'RoUi nde eTaI1l2006f) . 

I assume that MPHSs formed at redshift z ^ 3 and 
that they include LiBeB elements at the level of the 
IGM abundances at the time. If one can neglect the 
inhomogeneity of CCRs resulting from the local growth 
of a magnetic field, numbers of produced elements are 
proportional to target particle numbers. Consequently, 
resulting light clement abundance l/H docs not depend 
on density there. Although the metallicities like Fe/H 
of metal-poor stars reflect how the metal- enhanced and 
metal-deficient gas are mixed before star formation, the 
light element abundances are constant for the same for- 
mation epoch. In this case, there appear primordial 
plateaus on the plot of abundances as a function of metal- 
licity [Fe/H]. The energy fractions given to the CR ac- 
celeration, e = 0.31 for Model 1 and e = 0.029 for the 
rapid burst model realize the observed abundance of ^Li 
in MPHSs at z = 3. These fractions are reasonable, 
considering that the energy of Galactic CRs can be cov- 
ered by 10 — 30 percent of supernova remnant (SNR) 
energy (|Drurv et al.lll989f ). If '^Li observed in MPHSs 
is produced mainly by this CCR nucleosynthesis, these 
stars include Be and B which had been coproduced. 

In Fig. [4] the plateau levels of ^Li and ^Li in Model 
1 are drawn with observational data points in a plot 
of abundances as a function of metallicity [Fe/H]. As 
for the ^Li abundance, the BB N prediction i s cal- 
culated using the Kawano code (jKawanol Il992[ ) with 
the us e of the new world average of the neutron life- 
time (jMathews et al.l l2005[ ). I take the energy den- 



sity of baryons in the universe given by WMAP first 
year data analysis (jSpergel et al.ll2003h that is Vlbh'^ — 

0. 0224 ± 0.0009. This value corresponds to the baryon- 
to-photon ratio rj = 6.1 x 10~^°. The predicted abun- 
dance is then (''Li/H)BBN==4.5x 10~^°. The "^Li produc- 
tion by the CCR nucleosynthesis is a minor addition to 
the BBN result. The figure for the rapid burst model is 
very similar to Fig. |4l The factor of about three differ- 
ence between the BBN prediction and the observation of 
^Li abundance is apparent. 

The coproduced abundance levels of Be and B in Model 
1 are shown in Fig. [5] with observational data points. It is 
interesting that the predicted abundance levels of Be and 
B are located at nearly the lowest point ever detected. 
This model predicts primordial plateau abundances of 
Be and B as an analog of the likely ^Li plateau. If fu- 
ture observation of beryllium in metal-poor stars catches 
evidence of a plateau abundance, it would be explained 
by this CCRs origin. When the normal mode star for- 
mation alone is taken as the CR energy source, the re- 
sult is consistent with the observed abundance of ^Li 
in MPHSs at z = 3 with e — 0.73. This case corre- 
sponds to dashed lines in Fig. [5l In the same way, the 
coproduced abundances in the rapid burst model with 
e = 0.029 are shown as the dotted lines. This calcu- 
lation shows that Be and B abundances are somewhat 
higher in the rapid burst model than in Model 1 relative 
to that of ^Li. This calculation seems possibly to con- 
tain an error of the order of 10 % in Be and B abun- 
dances associated with the reading of a very brief burst 
of SFR (Fig. 13 of .Daigne e t al. 2006) and resulting sud- 
den e nrichment of C and O in the ISM of structures (Fig. 
16 of lDaigne et"aHl2006l . arXiv:astro-ph/0509183) as in- 
put profiles in this calculation. In the cosmic chemical 
evolution model of lDaigne et al.l (|200(l ). the rapid burst 
of star formation occurs in a very brief time and C and O 
abundances in structures increase accordingly refiecting 
the ejection of C and O nuclei by SNe. Since Be and B 
production is sensitive to the CR energy density times 
C and O abundances in structures, the fine-meshed time 
evolution of C and O abundances as well as CR energy 
density (which is associated with SN rate) is necessary to 
perform a precise calculation. However, since the rapid 
burst model gives sudden rises of C and O abundances, 
and C and O particles are accelerated before the ISM 
abundances can be somewhat diluted by accretion of the 
IGM, abundances of C and O in CCRs are relatively 
high. Then it is as expected that the rapid burst model 
produces more Be and B relative to ^Li than in Model 

1, considering that Be and B are produced by spallation 
processes of C and O while lithium isotopes are produced 
mainly hy a + a fusion. I show the results of calculations 
of light element production in the CCR nucleosynthesis 
model in Table [TJ In the second column, values of en- 
ergy fraction of SNRs given to CR acceleration required 
to produce ^Li at the observed level in MPHSs at z = 3 
are shown. The resulting abundances of light elements 
with the CR acceleration efficiencies in the second col- 
umn are listed in the third to seventh columns. 

In the calculations so far, the spectral index of the CR 
injection spectrum has been fixed to 7 = 3. I show the 
energy fraction of SNRs to CRs, e, required to produce 
^Li at the MPHSs level at z = 3 as a function of 7 in 
Fig. m The solid line corresponds to Model 1 and the 
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dashed hne to the case where only the normal mode stars 
are considered as the energy source of CCRs. One can 
see that the required energy fraction e is reasonable if 
2-7 7 ^ 3.1. Meanwhile, smaller spectral index of 
7 < 2.7 could not produce enough ^Li in this model. 

Figure [7] shows abundances of Be and B at z = 3 as 
a function of 7 in Model 1 with the value of e in Fig. 
[6l i.e., when ^Li is produced at the MPHSs level. It is 
found that relative produced abundances of Be and B are 
decreasing functions of 7. In other words, the steeper the 
momentum spectrum of CRs is, the more ^Li is produced 
relatively. In the region 2.7 < 7 < 3.1, where the ^Li 
abundance produced in this model can attain the MPHSs 
level with a reasonable partition of SNR energy to CR 

acceleration, log(Be/H) 13.1 - -12.9 and log(B/H)- 

— 12.0 — 11.8 are obtained. These abundances are again 
near the least abundances detected in metal-poor stars. 

Figure [8] shows the abundances of light elements pro- 
duced by the secondary process in the ISM as a function 
of redshift in Model 1 calculated by the second term in 
the rhs of Eq. (pij) . e = 0.31 is assumed to result in 
6Li/H=6 X 10-12 in ^j^g jq^ at 2 = 3. Solid lines cor- 
respond to the case where both the normal and massive 
modes are included, and dashed lines correspond to a 
contribution of the normal mode as the energy source. 
The contributions of primary and secondary processes to 
the light element synthesis in the ISM are estimated by 
the comparison between Fig. [1] and Fig. [51 Since ^Li and 
''Li are produced mainly by the a + a fusion which I con- 
tained in the calculation of Fig. [l] the contributions of 
the secondary process are relatively small. On the other 
hand. Be and B have some contribution from the sec- 
ondary process. For example, the Be and B abundances 
produced by the secondary process are about 25 % (65 %) 
of those by the primary process at z = 3 in Model 1 with 
both the normal and massive modes (with only the nor- 
mal mode). It is found that the BeB abundance levels 
produced by the OCR nucleosynthesis are roughly the 
same in the ISM and the IGM. This contribution from 
the secondary process is small in the rapid burst model 
(about 2.5 % for Be and B) at z = 3. The fraction of 
baryons included in the ISM is smaller at higher red- 
shift, when the energy injection in the rapid burst model 
occurs. It is then easily understood that the light ele- 
ment abundances produced by the secondary reactions 
are smaller in the rapid burst model considering the di- 
lution of light element abundances due to mass accretion 
by structure formation in the model. 

4. DISCUSSION 

The resulting abundances from Eq. (fT5)) have a clear 
dependence on input quantities as 

dNi }i{E, Z, Zs) , c / \rrCRi- ^I^IGM/ \ 

'—^^ dE oc 2_^£sn{zs)K'>p'^{zs)K^'^ (z). 

I adopted the results of the two models of Daign e et al.l 
llbOG) as the total CRs energy by SN explosions £sj^{zs) 
and the ratio of number abundance of i to that of p, i.e. 
i/p of CRs K^'{zs), which I assume to be the same as 
the ISM abundance of structures. These two quantities 
are uncertain, in fact. I use the primordial abundance for 
the background number ratio of j to p Kjp^{z). This 
quantity is reasonable and would not contain large uncer- 



tainty. A large difference in K^'{zs) {i=C,0) leads to 
a difference in produced Be and B abundances linearly, 
while ^Li and ''Li abundances have little influence since 
t hey are mainl y produced by the a + a fusion process. 

ISuzuki fc Yoshiil (j200lh have developed a model for 
the evolution of light elements in the Galaxy, which in- 
cludes the SN-induced chemical evolution with contribu- 
tions from SNe and CRs nucleosynthesis self-consistently. 
They have explained the light element abundances ob- 
served in metal-poor stars using their model with the 
CR abundances including SN ejecta 3.5 % in mass frac- 
tion. The linear relation between [BeB/H] and [Fe/H] 
has been obtained by the primary process to make Be 
and B. Consequently their CR abundance includes the 
C and O in mass fraction of (5 — 7) x 10^"^ originat- 
ing from the SN ejecta throughout the Galactic chemical 
evolution. On the other hand, the mass fraction of C-l-0 
is (3 - 10) X 10-3 in Model 1 of iDaigne et all (120061 ) 
I adopt which is roughly the same level of abundance 
as that ISuzuki fc Yoshiil ()2001[ ) used within a factor of 
< 2. The CR abundance I give in this study, there- 
fore, would be appropriate within a factor of 2, even 
if the real C and O abundances in the ISM producing 
CRs were lower than I give, supposing that the CRs have 
co ntribution from the SN ejecta by the fraction inferred 
bv ISuzuki fc Yoshiil (|200lh . 

I assume that all CRs escape from structures to the 
IGM, and do not consider a nonuniformity of the CR 
density in the universe. As structures grow in the uni- 
verse, a magnetic field grows accordingly, and the CR 
flux might get inhomogeneous, especially at low redshift, 
while the overproduction of ^Li pr ovides a constraint o n 
the confinement of CRs in the ISM (jRollinde et al.ll2006( l. 
Further study including the space distribution and time 
evolution of a magnetic field is desirable to estimate the 
light element abundances produced by the CCR nucle- 
osynthesis. 

The calculated result would also contain an uncertainty 
from two-step reactions, i.e., production of nuclide I2 by 
sequential non-thermal nuclear reactions: i + j — > h and 
h + 7 2 h with any nuc lid es h and .72 (iRamaty et al.l 
Il997t iKneller et all SoOSi) . IKneller et all (120031) have 



found that two-step reaction rates are only of the order 
of 1/10 smaller than one-step reaction rates. The effect 
of two-step reactions in the framework of this study is 
roughly estimated as follows. Nuclei produced by nu- 
clear reactions experience an energy loss and nuclear 
destruction. If one take nuclei with kinetic energy of 
£^ > 10 MeV/nucleon, the expansion loss is dominant 
loss process in the considered redshift range. The time 
scale for expansion loss is given by 



-^loss ' 



dE 



be^piE, z) 
E -f- En 



3 , n o M-1/2 



1-7 ~\~ rjc\ 

= 4.4xl0i's ° 



{E + 2E^) 



[an(i + z)3 + (i-a„)]" 



(23) 



where bc^p{E, z) is the energy loss rate for the cosmic ex- 
pansion. On the other hand, usi ng an empirical for mula 
of nuclear destruction (Eq. 1 in iLetaw et al.lll983t ). i.e.. 
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(T = 44. O^*''^ mb, the time scale for nuclear destruction 
is given by 

= 8.0 X 10^0 s f—j + (24) 

Ratios of time scales for nuclides with ^ ~ 10 and E — 
10 - 10^ MeV/nucleon are Tioss/Tb 10"^ - 10-2 at 
z = 3, 3 X 10-3 - 3 X 10-2 at z = 10, and 10"^ _ o.2 
at z = 30, respectively. Therefore, the fraction at which 
nuclei produced by nuclear spallations experience second 
nuclear reactions before losing enough energy is expected 
to be at most of the order of 10 %, if the energy 
spectrum of CRs indicates small amount of high energy 
CRs. I check a fraction of high energy CRs in number, 
which would be estimated by 



100 MeV 



Q,{E',z)dE' _ P(100 MeV) 



LfMlvQ'^iE',-)dE' P(lOMeV) 



(25) 



where I defined 



P{E) 



E' + Eq 



E [E'{E' + 2Ea' 



l(7+l)/' 



^^dE'. (26) 



For 7 = 2,3,4, F(2) = 0.31, F{3) = 0.095, and i^(4) = 
0.029 are obtained, respectively. The fraction of high 
energy CRs is thus relatively low for the range of 2 < 
7 < 4. As a result, the effect of two-step reactions is 
small and would be at most of the order of ~ 10 %. 

.RoUi nd e et al.l jlOOS) also study the CR production 
of Be and B by CCRs. Their conclusion is very sim- 
ilar to that of this study, and a potentially detectable 
Be and B is produced by CCR-induced spallation reac- 
tions at the time of the formation of the Galaxy (z ~ 3) . 
However, there are some differences of assumptions be- 
tween the two studies, which are compared here. I give 
the CR spectrum of CO nuclides by the same shape as 
those of p and a particles, while iRollinde et al.1 ()2008D 
give it by a broken power law which matches the ob- 
served present-day Galactic CR spectrum. Moreover, 
I give the abundances of CR by those of structures in 
the co smic chemi cal evolut ion model of iDaigne et al.l 
(|2006f l. while RoUi nde et all 



2008fl give them by those 



of structures in .Daigne et al.l ( 20061 ) model multiplied by 
abundance enhancement factors of present Galactic CR 
fluxes. I assume that the CR confinement by a mag- 
netic field is ineffective in the early universe, and that 
all CRs ge nerated by SNe in str uctures escape to the 
IGM, while iRolhnde et"all ()2008D apply a shape of CR 
diffusion coefficient, which leads to some degrees of the 
nuclear destruction of CRs in structures and the applica- 
tion of a broken power law in the CR energy spectrum of 
CO nuclides. Their diffusion coefficient is derived assum- 
ing that it is given by that in a magneto-hydrodynamic 
(MHD) turbulence, and that the magnetic energy density 
is proportional to the thermal energy, and that a char- 
acteristic length scale for the magnetic field is given by 
the local Jeans scale. It is interesting that the two stud- 
ies using different assumptions for the uncertain physical 
inputs arrived at the similar conclusions. 

I comment on a difference between the results of this 
CCR production of the l ight elements and those o f 
the flare production model (jTatischeff fc Thibaiidll2007[ ) . 



The CCR nucleosynthesis leads to the production of Be 
and B at the lowest level ever detected, in this model 
calculation. The nucleosynthesis on the main sequence 
stars triggered by flare-accelerated nuclides, on the other 
hand, results in the negligible productions of ^Li, Be and 
B, which are proportional to the metallicity and exist 
at very low abundance levels under the observed linear 
trend as a function of metallicity. Therefore if we observe 
a signature of the primordial origin of Be and B by mea- 
surements of MPHSs, the production mechanism of Be 
and B would not be the flare-energized nuclear reactions 
after the star formations, and it would be thought that 
the plateau abundances of Be and B originate in the CCR 
production and ^Li has been coproduced at pregalactic 
phase by the CCRs. 

5. CONCLUSIONS 

The recent observations of MPHSs reveal the proba- 
ble existence of high plateau abundance of ^Li, which 
is about a thousand times higher than predicted in the 
standard BBN model. Since the standard Galactic chem- 
ical evolution model with the Galactic CR nucleosynthe- 
sis gives lower values of ®Li abundance at the metallici- 
ties of the observed MPHSs, some mechanism must have 
produced ^Li existing in the surface of MPHSs. As a can- 
didate of a ^Li productio n mechanism, the ear ly burst of 
CRs has been proposed (RoUin de et al.l 120051 ) . and the 
nucleo synthesis by the CRs from SN explosions is calcu- 
lated (|RoUinde et al.ll2006D in a detailed model of cos- 
mic chemical evolution llPa igne et al.ll2006f ) which satis- 
fies various observational co nstraints includ i ng an early 
reionization of the universe. IRollinde et all (|2006D have 
found that the a + a fusion reaction can produce ^Li to 
the level observed in MPHSs. 

I calculate the cosmological cosmic ray nucleosynthesis 
of Be and B isotopes as well as ^Li and ''Li with the use 
of Mo del 1 and the rapid burst model in iDaigne et al.l 
(j2006l ). It is assumed that all CRs produced by SNe in 
the ISM escape to the IGM and the CR intensity is al- 
ways homogeneous in the universe. I found tha t when 
Model 1 (the rapid burst model) of IDaigne erall (|2006l ) 
is adopted for the SFR of the universe and the metal 
abundances of CRs, Be and B are produced at (above) 
the levels observed in the most metal-poor stars with 
detection of Be or B, if the ^Li plateau abundance is 
made by the same CCR nucleosynthesis. The CR ac- 
celeration energy needed to make ^Li primordial plateau 
abundance at the observed level is ~ 3 — 31 % of the 
SN kinetic energy. This value is not too large in view 
of an inferred present fraction of the SN energy used to 
the CR acceleration (|Drurv et al.lll989l ). The pregalac- 
tic SN activity might have produced some level of light 
elements. Although the resulting abundances of light el- 
ements d epend on the param eters which I fixed to the 
values of IDaigne et al.l (j2006f ). the future measurements 
of metal-poor stars would show the reasonableness of this 
early ^LiBeB production mechanism, and might provide 
a signature of a primordial Be (and perhaps B) plateau. 
Further observations of LiBeB elements in MPHSs are 
highly desirable and valuable. 

I thank Frederic Daigne for explaining the result of his 
chemical evolution calculation, and Emmanuel RoUinde 
for helpful comments. I am grateful to Adrian Jenkins 
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Fig. 1. — Abundances of light elements in the IGM as a function of redshift in Model 1 (solid lines), e = 0.31 is assumed to result in 
^Li/H=6 X 10~^^ at z = 3. The contribution of the normal mode stars only to the light element production is shown by the dashed lines. 
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Fig. 2. — Abundances of light elements in the IGM as a function of redshift in the rapid burst model, e = 0.029 is assumed to result in 
8Li/H=6 X 10-12 at z = 3. 
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Fig. 4. — Plateau abund ances of lithium isot opes produced b y the OCR nucl e osynt hesis in Model 1 wi th the accelerat ing efficiency 
e = 0.31. '^Li data are fromlAsplund et al.l II2006I . filled triangles'). IBonifacio et al.l 1)20071 . open squares), and lShi et al.l l|2007l . open stars). 
^Li data are from lAsplund ct al. (200^, large filled circles to detections, small filled circles to upper limits) and llnoue "etaTI ||2005| . open 
circles to detections, crosses to upper limits). 
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Fig. 5. — Plateau abundances of Be and B produced by the CCR nucleosynthesis in Model 1 (solid lines) with the accelerating efficiency 
e = 0.31. The case that the normal mode star formation alone is considered as the CR energy source corresponds to the dashed lines 
with e = 0.73 to realize the MPHS value of ®Li at z = 3. Plateau abun dances in the rapid burst m odel with the accelerating efficiency 
e = 0.029 are also sh own as the dotted li nes. ^Be data are from Bocsgaar d et al] I I1999I . filled circles) . IPrimas et al. (2000a, open circles), 
Primas et all Il2000b|. ope n triangle), and Boesgaard & Novicki (2006, open stars) . B data are from [ Duncan et al.l 11199 TL filled squares), 
G^ci^^opez~et''an " 1998 . open squares). .Primas et al.. (1999. . filled triangles), and lCunha et ahl II200C , crosses). 
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Fig. 6. — Energy fraction of SNRs to CRs, e, as a function of the index of the CR injection spectrum, 7, required to produce ®Li at the 
MPHSs level ^Li/H=6 X 10~^^ at z = 3. The solid line corresponds to Model 1 and the dashed line to the case where only the normal 
mode stars are considered as the energy source of CCRs. 
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Fig. 8. — Abundances of light elements produced by the secondary process in the ISM as a function of redshift in Model 1 (solid lines), 
e = 0.31 is assumed to result in ^Li/H=6 X 10~^^ in the IGM at 2: = 3. The contribution of the normal mode stars only to the light element 
production is shown by the dashed lines. 
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TABLE 1 

Abundance results for two models 



Model 


e 


^Li/H 




''Li/H 




^Be/H 




i«B/H 




"B/H 




Model 1 


0.31 


6.0X 10" 


-12 


7.3x 10" 


-12 


9.0X 10" 


-14 


3.2x 10" 


-13 


7.6x 10" 


-13 


Model 1 (Pop II only) 


0.73 


6.0X 10- 


-12 


7.3x 10- 


-12 


5.7x 10- 


-14 


2.1x 10- 


-13 


4.9x 10- 


-13 


Rapid burst model 


0.029 


6.0X 10" 


-12 


7.9x 10" 


-12 


4.4x 10" 


-13 


2.0x 10" 


-12 


5.2x 10" 


-12 



